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Abstract 


Conser'  ation  of  potential  vorticity  in  Eulerian  fluids  reflects  particle  interchange 
symmetry  in  the  Lagrangian  fluid  version  of  the  same  theory.  The  algebra  associated 
with  this  symmetry  in  the  shallow  water  equations  is  studied  here,  and  we  give  a 
method  for  truncating  the  degrees  of  freedom  of  the  theory  which  preserves  a  maxi¬ 
mal  number  of  invariants  associated  with  this  algebra.  The  symmetry  associated  with 
keeping  N  modes  of  the  shallow  water  flow  is  SU(N).  In  the  limit  where  the  number  of 
modes  goes  to  infinity  {N  —*  oo),  all  the  conservation  laws  connected  with  potential 
vorticity  conservation  are  recovered.  We  also  present  a  Hamiltonian  which  is  invari¬ 
ant  under  this  truncated  symmetry  and  which  reduces  to  the  familiar  shallow  water 
Hamiltonian  when  N  — ►  oo.  All  this  provides  a  fiiiite  dimensional  framework  for  nu¬ 
merical  work  with  the  shallow  water  equations  which  preserves  not  only  energy  and 
enstrophy  but  all  other  known  conserved  quantities  consistent  with  the  finite  number 
of  degrees  of  freedom.  The  extension  of  these  ideas  to  other  nearly  two  dimensional 
flows  is  discussed. 
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1  Introduction 


Many  geophysical  problems  are  naturally  decomposed  into  a  many  layered  approximation 
with  each  layer  governed  by  the  so-called  shallow  water  equations  [l].  These  equations  take 
the  fluid  density  to  be  constant  in  each  layer,  and  because  the  horizontal  dimensions  are 
assumed  much  larger  than  the  vertical,  hydrostatic  balance  is  taken  to  hold  in  each  layer 
separately.  Vertical  variations  in  each  layer  are  ignored  in  the  two  dimensional  horizontal 
velocity  v(x,<)  =  [t;i{x,y,t),t;2(a:,y,t)]  (x  =  (i,y)),  and  incompressibility 
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determines  the  vertical  velocity  ^^(x  zU).  Using  local  hydrostatic  balance,  the  pressure  is 
eliminated  in  terms  of  the  thickness  of  the  vertical  layer,  h{x,t),  and  h{x,t)  becomes  the 
third  dependent  dynamical  variable  for  the  reduced  system. 

The  evolution  equations  for  v(x,t)  and  /i(x,t)  serve  both  as  a  useful  model  for  the  dy¬ 
namics  in  a  thin  layer  of  fluid  and  as  an  important  ingredient  in  more  complicated  models  of 
the  whole  atmosphere  or  ocean  [2j.  In  complex  models  which  attempt  to  represent  the  full 
dynamics  of  the  atmosphere,  for  example,  one  must  add  to  the  basic  shallow  water  equa¬ 
tions  additional  dynamics  to  describe  radiative  transfer,  internal  waves,  cloud  formation, 
relevant  chemistry,  etc.  Whatever  the  goal  of  the  dynamics  of  the  shallow  water  equations, 
if  one  is  to  solve  these  equations  some  form  of  truncation  of  the  infinite  degrees  of  freedom 
must  be  made  to  progress  numerically.  Truncations  directly  in  Eulerian  or  Lagrangian 
configuration  space  or  in  the  dual  Fourier  space  fail  to  preserve  all  the  conservation  laws 
respected  by  the  underlying  particle  interchange  symmetry  of  the  Lagrangian  theory  which 
exhibits  itself  in  the  conser'^^tion  of  potential  vorticity  in  Eulerian  fluid  mechanics.  These 
latter  remarks,  of  course,  apply  only  when  the  shallow  water  flow  is  inviscid,  as  we  shall 
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assume  throughout  this  paper.  We  shall  have  a  few  remarks  to  make  at  the  end  of  this 
paper  about  the  potential  use  of  our  results  for  the  case  with  friction. 

In  this  paper  we  take  up  the  much  studied  subject  of  shallow  water  equations  with  the 
goal  of  providing  a  truncation  of  the  degrees  of  freedom  from  infinity  to  a  finite  number 
using  a  method  which  preserves  the  maximum  number  of  conserved  quantities  consistent 
with  the  reduction  in  the  number  of  degrees  of  freedom.  When  this  number  returns  to 
infinity,  that  is  when  the  truncation  is  removed,  the  theory  preserves  all  the  quantities 
associated  with  potential  vorticity  conservation.  Our  work  takes  place  in  the  Lagrangian 
formulation  of  the  theory.  The  truncation  is  made  in  the  Fourier  space  of  variables  dual 
to  the  Lagrangian  labels  of  fluid  particles.  The  algebra  associated  with  the  symmetry  of 
pai  tide  interchange  is  altered  as  part  of  the  truncation,  and  the  finite  number  of  Casimir 
invariants  of  the  new  algebra,  which  is  SU(N)  and  thus  familiar,  replace  the  infinite  number 
of  conserved  quantities  following  from  potential  vorticity  conservation.  In  the  limit  N  — > 
oo,  the  usual  conserved  quantities  are  recovered. 

The  methods  we  shall  use  we  first  learned  from  two  sources.  One  is  the  work  by  Fairlie, 
et  al  [3]  on  finite  algebras  in  string  theory,  and  the  other  is  an  application  of  those  methods 
to  the  two  dimensional  Euler  equations  independently  invented  by  Rouhi  [4]  and  Zeitlin  [5]. 
The  latter  application  may  be  quite  interesting  in  other  geophysical  applications  where  two 
dimensional  Eulerian  flows  are  studied,  but  we  have  not  pursued  that  line  of  investigation. 
We  have  analyzed  the  shallow  water  equations,  as  presented  here,  both  for  their  interest 
as  indicated,  and  also  because  they  have  numerous  useful  formal  similarities  with  internal 
wave  dynamics  and  with  surface  wave  physics.  We  will  indicate  these  similarities  at  the 
end  of  this  paper.  Our  work  here  is  also  in  planar  geometries.  The  extension  to  flows  on  a 
sphere,  while  algebraically  complicated,  is  more  or  less  straightforward  in  concept  as  seen 
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in  the  paper  of  Hoppe  [6]. 

The  shallow  water  equations  and  their  numerical  solution  using  various  truncations  has 
become  a  subject  of  renewed  interest  of  late  because  of  the  work  to  place  these  equations 
and  their  more  complex  forms  on  parallel  processing  machines  (7j.  The  goal  of  that  ef¬ 
fort  is  to  build  numerically  efficient  climate  models  for  investigations  of  very  long  times 
(thousands  of  simulated  years)  and/or  issues  requiring  very  high  spatial  resolution.  We 
expect  that  the  truncation  presented  here,  which  by  its  formulation  preserves  as  much 
of  the  original  symmetry  as  possible  of  the  structure  of  the  shallow  water  equations,  will 
prove  an  attractive  alternative  to  straightforward  finite  element,  discrete  spatial  grid,  or 
spectral  methods  for  these  equations. 

In  the  next  Section  we  review  the  shallow  water  equations  in  Eulerian  and  Lagrangian 
formulation  and  write  down  the  algebraic  structure  associated  with  particle  interchange 
symmetry.  The  third  Section  is  devoted  to  the  SU(N)  truncation  of  the  theory  in  La¬ 
grangian  formulation  and  also  presents  the  truncated  Hamiltonian  for  the  shallow  water 
flow.  The  final  Section  has  our  comments  about  further  uses  of  our  observations  in  otlier 
problems  of  geophysical  interest  and  contains  the  summary  of  our  present  work. 

2  Shallow  Water  Theory 

2.1  Equations  of  Motion  and  Symmetry 

The  Eulerian  shallow  water  equations  govern  the  evolution  of  a  two  dimensional  horizontal 
velocity  v(x,t)  in  a  fluid  of  vertical  thickness  h(x,t}  of  the  fluid  via  the  familiar  evolution 
equations 

+v(i,t)  •  Vv(x,t)  =  -gVh(x,t) 
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+ V-[/i(x,y)v(x,01  =  0, 


(2) 


where  g  is  the  gravitational  constant  and  V  =  \di,d-i\.  If  the  frame  is  rotating  about  the 
z-axis  at  angular  velocity  //2,  a  term  v{x,t)  x  i/(x)  appears  in  the  equation  for  v{x,i). 
As  indicated  above  these  equations  follow  from  the  three  dimensional  Euler  equations  of 
a  thin,  homogeneous  fluid  with  hydrostatic  balance  in  each  layer  determining  the  pressure 
p{x,zj)  in  that  layer  in  terms  of  the  thickness  p{x,^,^)  =  g[h{x,t)  —  z\. 

The  total  energy 

W£(v,/.)  =  i/  i=x||v(x,l)r+Mx,fn,  (3) 


is  conserved  by  solutions  to  these  equations  and  the  Eulerian  potential  vorticity 


z  •  V  X  v(i,  t) 
h{x,t) 


(4) 


satisfies 

^5^^  +  v(x,()-g£(x,()=0,  (5) 

so 

j  d^xh(x,t)G(q£(x,t))  (6) 

is  time  independent  for  arbitrary  G(g£). 

These  conservation  laws  arise  from  the  particle  interchange  symmetry  exhibited  by 
the  canonical  or  Lagrangian  formulation  of  the  theory.  In  Lagrangian  formulation,  the 
dynamical  variables  are  the  particle  position  Y(r,t)  at  every  particle  label  r  and  time  and 
the  conjugate  momentum  P(r,t).  In  terms  of  these  variables  the  Hamiltonian  reads 


Wi(Y,P)  =  i/  <i“rl|P(r,l)l“+ff^(Y(r,())-'|,  (7) 
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where  the  Jacobian 


J(Y(r,0)  = 


a(Y) 


has  been  introduced.  In  these  equations  and  below  a  Lagrangian  density  po(r)  wliich 
describes  the  distribution  of  mass  among  the  particles  labeled  by  r  has  been  absorbed  into 
the  definition  of  the  labels  themselves  [8). 

The  evolution  in  time  of  any  functional  yl[Y(r,t),P(r,t)]  follows  from  Hamiltonian’s 


equations 


=  {^|Y(r,O.P(r,l)l,/f(Y,P)}. 


introducing  the  canonical  Poisson  brackets  between  functionals  yi(Y,P)  and  5(Y,P) 

The  fundamental  canonical  brackets  between  Y(r,t)  and  P(r,t) 

-t')  ( 


follow  from  this  bracket,  and  we  shall  use  it  extensively  below. 


In  this  notation  the  potential  vorticity  takes  the  form 


nfr-  t\  -  f  9yaiT^,i)dPa{r,t) 
9(^1  V  ^ab  ^  n  ’ 


where  the  indices  a, 6, a,. ..  run  over  1  and  2,  and  repeated  indices  are  summed  over,  tat 
is  the  completely  antisymmetric  symbol  in  two  dimensions.  The  conservation  law  is 


and  all  derivatives  in  Lagrangian  formulation  are  taken  at  fixed  r.  The  integrals 


C|Y,P|=|  <iVG(,(r.()) 
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are  clearly  constant  in  time  for  any  <^(9). 

The  generators  of  the  particle  interchange  symmetry  are  determined  by  examining  the 
invariance  of  the  action  defining  Hamilton’s  principle 

ACTION  =  £■  /  d‘r  -  SJ-I,  (15) 


under  variations  in  r  at  fixed  Y  (r,  t)  and  fixed  J.  These  variations  must  be  of  the  fo»-m  (9] 

dA{T) 


SVa  —  Ca6" 


drh 


(16) 


with  >l(r)  an  arbitrary  function  of  r.  The  time  independent  functional  C[Y,P)  defined 
above  has  the  property  that 

=  (17) 

which  is  of  the  correct  form.  So  the  global  generators  of  the  particle  interchange  symmetry 
are  just  integrals  of  arbitrary  functions  of  the  potential  vorticity.  The  algebra  associated 
with  the  particle  interchange  symmetry  is  rather  rich.  The  potential  vorticity  does  not 
Poisson  commute  with  itself,  but  the  Poisson  brackets  of  the  9{r,t)  do  close  on  themselves. 
We  turn  to  this  algebra  now. 

2.2  Potential  Vorticity  Algebra 

To  study  the  algebra  of  the  potential  vorticity  we  have  found  it  most  straightforward  to 
work  in  the  Fourier  space  dual  to  the  particle  labels  r.  We  continue  to  work  in  planar 
geometry  and  limit  horizontal  space  to  a  box  of  size  L  x  L,  In  this  box  we  introduce 
Fourier  transforms  as  follows:  for  functions  /(r)  we  write 


00 

/W  =  <7(n)exp{{Kn-r|, 

n=— 00 


(18) 
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and  the  inverse 


s(n)  =  /W  exp[-i/cn  •  r],  (19) 

where  the  vectors  n,m, . . .  are  composed  of  integers  n  =  [nijUj];  ria  =  0,  ±1,±2, . . and 
K  =  21^ I L.  We  take  the  Fourier  decomposition  of  the  canonical  variables  to  be 

=  7ZlQ(n)exp[i/cn-r] 

^  n 

P(r,0  =  ~Y^n{n)exp[{Kn-T],  (20) 

^  n 

which  gives  the  Poisson  brackets  in  Fourier  space 

{Qa(ii),  nj,(m)}  =  (21) 

Using  this  Fourier  transformed  set  of  canonical  variables  we  have  for  the  potential  vorticity 

9(r)  =  ^^X^?(n)exp[i'/cn  t),  (22) 

^  n 

with 

m 

We  have  defined  the  quantity  nxm  =  nxmz  —  n^mi  is  this  equation,  and  the  normalization 
for  $^(n)  has  been  chosen  to  make  the  last  formula  and  many  to  follow  ais  simple  as  possible. 

The  algebraic  properties  of  f  (n)  arc  seen  in  the  Poisson  brackets 

{?(n),Qa{ni)}  =  nxmQa{n  +  m), 

{f(n),na(m)}  =  nxmna(n  +  m),  (24) 

and 


{f (n).  d™)}  =  n  X  m  f (n  +  m). 


(25) 


The  last  Poisson  bracket  exhibits  the  structure  of  the  particle  interchange  algebra  and 
shows  it  to  be  an  infinite  dimensional  algebra,  as  we  might  expect  since  it  comes  from  a 
continuum  set  of  operations  on  label  space  r.  This  set  of  continuum  operations  translates 
via  the  Fourier  transform  into  a  discrete  infinity  of  operation  in  n  space. 

We  will  also  require  the  properties  of  the  Jacobian  under  Poisson  commutation  with  the 
potential  vorticity  f(n)  since  we  will  be  interested  in  construction  a  Hamiltonian  invariant 
under  the  particle  interchange  symmetry.  The  Hamiltonian  of  the  shallow  water  equations 
is  invariant  under  this  symmetry,  which  is  our  starting  point  for  all  this  after  all,  and  this 
invariance  is  the  same  statement  as  the  time  independence  of 

^^  =  {?(r.<),ff(Y,P)}=0.  (26) 

When  we  truncate  the  Fourier  transform  variables  in  just  a  moment  we  will  want  to  retain 
this  invariance  under  the  truncated  potential  vorticity.  The  transformation  properties  of 
the  Jacobian  under  f(n)  will  be  of  interest  in  this. 

The  Fourier  decomposition  of  the  Jacobian  is 

r,i))  = 

=  ^^IIp(n)exp(f«n-r],  (27) 

^  n 

which  gives 

pM  =  ^  X  mQ(m)  X  Q(m')«5n.m+m' 

m,m' 

=  X  m<5i(m)Q2(n  -  m).  (28) 

m 

The  Poisson  bracket  of  f(n)  and  p(m)  is  then 

{f(n),p(m)}  =  n  X  mp(m-f-n).  (29) 


djY) 

d(r) 
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One  may  summarize  these  brackets  as  saying  that  under  the  infinite  algebra  of  the 
potential  vorticity  f(n)  the  anaraities  Q(m), /T(m),  jo(rn)  and  c{rn)  itself  transform  as 
“vectors”. 

3  TriTncation  of  the  Modes 

3.1  Algebra  of  the  Truncated  Modes;  SU(N) 

Our  truncation  of  the  particle  interchange  algebra  is  motivated  by  the  idea  of  limiting 
the  number  of  Fourier  modes,  but  a  simple  cutoff  on  the  components  of  the  vectors  u 
does  not  respect  the  Poisson  brackets.  Suppose  we  limit  each  component  of  our  integer 
vectors  n  =  [ni,n2]  to  —A/  <  <  M,  then  the  presence  of  the  sum  vectors  n  +  ni 

in  the  Poisson  brackets  of  f(n)  means  that  vectors  in  the  range  are  mapped  out  of  the 
range.  To  address  this  we  note  that  restricting  vectors  in  label  space  to  a  box  of  size  L  x  L 
and  making  the  Fourier  transform,  we  have  in  effect  mapped  our  space  onto  a  torus  by 
implicitly  identifying  the  sides  under  the  assumption  of  periodicity.  If  we  were  to  formalize 
this  periodicity  by  requiring  all  sums  of  integer  vectors  to  lie  within  the  range  [  — Af,  Afi  by 
a  modulo  or  remainder  operation,  we  would  still  need  to  deal  with  the  terms  n  x  ni  that 
appear  in  all  of  the  Poisson  brackets  with  These  are  the  so  called  structure  constants 
associated  with  the  group  properties  of  particle  relabeling  invariance,  so  it  is  suggestive 
that  modifying  them  cis  well  would  be  required  to  make  a  consistent  theory  of  truncated 
Fourier  modes.  In  particular  one  must  address  the  Jacobi  identity  which  provides  the 
statement  that  the  operations  in  question  do  close  to  form  a  group. 

The  set  of  operations  which  provides  a  consistent  truncation  of  the  modes  comes  from 
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changing  the  definition  of  the  Fourier  components  of  9(r,t)  by 

fyv(n)  =  “  Yl  sin[/C7»,n  x  mjna(m)Qa(n  -  rn),  (30) 

IT1=-Af 

where  N  =  2M  +  1,  all  components  na,mat...  are  restricted  to  [-M,M],  and  =  2njN. 
In  the  limit  A/  — >  oo  or  equivalently  N  —*  oo,  namely  as  /c//  — »■  0,  this  definition  of  the 
potential  vorticity  is  the  same  as  in  the  original  Fourier  transform.  The  definition  of  the 
Fourier  density  for  the  Jacobian  is  modified  to 

1  M  1 

PN{n)  =  ~  Y  —  sin{KAfm' X  m)Q(m)  X  Q(m')(5n.m+m' 

^  xn.Tn'=~M 

M  j 

=  Y  — sin[x/vn  X  m|  Qi(m)Q2(n  —  m).  (31) 

m=-Ar 

The  definitions  of  the  Q(n)  and  P(n)  are  unchanged  and  the  Poisson  brackets  between 
them  are  still 

{Fo(n),  n(,(m)}  =  <5a(,(5o,m+n>  (32) 


with  the  rule  that  vector  components  out  of  (-M,M]  are  mapped  back  into  the  range. 

Now  the  Poisson  brackets  among  the  other  variables  Y(n), /7(rij,  pAf(n) 

are  found  to  be 


Y  AT 

{<rN(n)i<?a(m)}  = - Y  sinf/c/^n  X  m'|(5o,ni'+mQo(n  -  m') 

m’=-Af 

=  —  sin[/c;^n  X  ni]i5a(n  +  m), 

kn 


(33) 


and 


{?N(n),na(m)}  = 

{?N(n).P/v(m)}  = 

{fiv(n),fiv(m)}  = 


—  sin[/c/fn  x  mjno(n  +  m), 

—  sinfKArn  x  m]/9Ar(n  +  m), 
Kn 

—  sin[/c;vn  x  ni). 


(34) 

(35) 

(36) 
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In  each  of  these  Poisson  bracket  relations  we  have  used  the  trigonometric  identity 
sin(a[m  x  m'  —  m  x  n])  sin(a[n  x  m'|)  -f 

sin(a[n  x  m  —  n  x  m'j)  =  sin(a(n  x  m])  sin(a({n  +  m)  x  m'l),  (37) 
for  an  arbitrary  (complex)  constant  a.  For  us  a  — 

This  set  of  Poisson  brackets  defines  a  algebra  which  is  SU(N)  with  the  f(n)  as  generators 
of  infinitesimal  SU(N)  transformations.  The  truncation  of  the  Fourier  modes  with  the 
modification  of  the  Poisson  brackets  now  provides  a  consistent  reduction  from  an  infinite 
number  of  modes  to  N.  The  critical  issue  in  checking  this  consistency  is  verifying  that 
the  Jacobi  identity  among  Poisson  brackets  is  satisfied,  and  with  the  change  of  structure 
constants  from  n  x  m  ~+  ^sin[/c;y^n  x  m]  this  is  readily  established. 

We  want  to  construct  a  Hamiltonian,  Hs,  now  which  is  invariant  under  this  SU(N) 
and  becomes  the  shallow  water  Hamiltonian  in  the  limit  N  oo.  For  this  purpose  we 
need  to  investigate  the  Casimir  invariants  o^  the  SU(N)  algebra.  These  we  generalize  from 
the  observations  of  Fairlie  et.  al.  [lO].  If  the  Poisson  brackets  of  ?jv(n)  with  any  function 
in  Fourier  space  /(ni)  is 

{fw(n)./(in)}  =  —  sin{/c;vn  x  mj/(n  +  m),  (38) 

Kff 

then  the  following  sums  Poisson  commute  with  f/v(n) 

M 

=  Z]  /(n)/{m)^o.n+m, 
n,m=-M 

M 

exp[z'«:/v(ni  x  nj  +  nj  x  ns  +  nj  x  n3j<5o,n,+n3ni;^9) 

n,--M 

and  similarly 

M 

Ci+iC/)  =  IT  n  exp[i/c/vnc  x  n^|/(n,)/{n2) . . .  /(nt)/(-(ni  +  ns  +  . . .  +  n^)) 

ni,...,nx,=  *4  a<.0 

(40) 
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These  sums  go  over  in  the  limit  — ♦  0  to  the  integrals 


/  (41) 

for  y  =  1, . . .  L.  We  will  use  this  in  constructing  Hs- 

3.2  The  Truncated  Hamiltonian:  Hff 


We  begin  by  remembering  how  the  full  Lagrangian  variable  Hamiltonian 

H[Q{n),n(n))  =  \  £  |n,(n)n.(-n)|  +  9  / dV  J’",  (42) 

^  n--co  •  ■' 

conserves  f(n).  The  Poisson  brackets  of  f(n)  with  the  first  term  is 
1  “ 

-  [n  X  mna(m  +  n)na(-m)  -  n  X  mna{m)na(n  -  m)],  (43) 

^  m=-oo 

which  vanishes  by  symmetry;  let  m  —*  -m  in  the  second  sum.  The  Poisson  bracket  of  the 
second  term  is 


.2  5{exp[-i«;n  •  r],  J~^) 
o  r  — 


5(n,r2) 


(44) 


which  vanishes  upon  integration  by  parts. 

When  we  change  the  structure  constants  in  the  Poisson  brackets  and  truncate  the  modes 

so  all  n  are  in  [-M,M],  the  Poisson  bracket  of  fjv(n)  with  the  kinetic  energy  becomes 

1  M  M 

£  na(m)na(-m)}  =  —  ^ 

^  m=-M  m=-Af 


(sin[/c;vn  x  m]/7(n  +  m)  •  /7{— m) 

—  sin[K^/n  X  m]JrZ'(m)  •  /7{n  -  m)],  (45) 


which  again  vanishes  because  of  symmetry. 
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The  Jacobian  term  is  more  tricky.  The  Poisson  bracket  of  with 


P/tf(m)exp[iKyvn-m|, 


(46) 


n=-M 


IS 


(f  (n),  J;y(m)}  -  exp[-i«;vn  ■  m](J;^(m  +  2  x  n)  -  -  2  x  n)).  (47) 

4m 

Though  in  the  limit  as  k//  — ♦  0,  this  becomes  the  same  as  the  continuum  theory,  it  shows 
that  a  Hamiltonian  constructed  out  of  integrals  over  tf^r  of  will  not  Poisson  commute 
with  f{n). 

We  choose  instead  to  represent  the  quantity  J”*  which  appears  under  the  integral  in 
the  Hamiltonian  as  an  expansion  in  powers  around  J  =  Jq.  The  full  expansion  is 

1 


1 

J 


Jo  -  {Jo  ~  J) 

1  7 


which  we  truncate  at  /c  =  TV  —  1  to 


'0  k-O 


1  1  J 

7  -  IE(i-t) 


J 


(48) 


^0  t=o 


Jo 


(49) 


Clearly  for  (l  —  X  this  expression  is  essentially  j .  It  is  also  a  hnite  sum  of  powers 

of  J,  and  this  leads  to  the  approximation  we  will  use  in  H/^: 

1 


7 


1  ^  J 

=  TE(-f)‘"c^ 

Jo  Jo 


(50) 


which  should  be  good  for  |1  ~  ;^|  <  1  or  0  <  ;^  <  2,  which  will  suit  our  purposes  quite 


well. 
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The  idea  is  then  to  approximate  /  <i*r  j  by 


jtc;  (51) 

Jq  J  Jq 

and  then  to  replace  ^  d'^r  by  times 

M 

Cp+i{PAr)  =  £  n 

ni....,np=-M  a</S 

exp(i/Civna  x  n;j]pAf(ni)pA,{n2) . .  .pA/{np)pjv(-(ni  +  nj  +  . . .  +  np)),(52) 

recalling  that 

J'N{n)expl-ix/vm-n].  (53) 

(ZTTJVJ  jj=-A/ 

The  Hamiltonian  invariant  under  the  action  of  the  local  SU{N)  generators  f  (n))  is  then 

1  M  N  (2Tr\2k-2 

Hx  =  5 n{^)  •  ^r{-n)  +  ^  E (54) 
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